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Abstract 

^ , Transition from Fourier series to Fourier integrals is considered and 

error introduced by ordinary substitution of integration for summing is 
estimated. Ambiguity caused by transition from discrete function to con- 
tinuous one is examined and conditions under which this ambiguity does 
not arise are suggested. 



^ ■ 1 Introduction 

-4— > 

Quite often it appears to be appropriate to solve some physical problems at 
pi il first in a finite size box (see e.g. 012113) and to increase the volume infinitely 

■ only at the final stage. In particular, such approach provides a reasonable 

regularization. As a rule periodic border conditions are imposed, so Fourier 
transform in spatial variables becomes a discrete function (Fourier coefficients), 
which is assumed to convert into a continuous one, when borders are pulled 
apart and period 2t:t is infinitely increased. Ordinary substitution integration 
' for summing is acceptable for 'smooth enough' and 'fast enough' decreasing 

functions, but being applied to distributions it may introduce errors which we 
estimate in this paper. Moreover, a transition from a discrete function to a 
continuous one may introduce ambiguity. We estimate possible ambiguity and 
consider conditions under which such ambiguity doesn't arise. It should be 
noted that even if the original function is smooth, after transition r — > cxd it 
may convert into a distribution. 

In this paper we confine ourselves to the case of tempered distributions 
F {x) e S", which are defined as functionals {F (x) , cj) {x)) on fast decreasing 
test functions (j>{x) G S ^ElEl- Recall that F (x) is a tempered distribution, if 
and only if it is a finite order derivative of some continuous tempered function 

*E-mail address: vkpotrov@yandex.ru 
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G (x), i.e. F [x) = G*^"^ (x), n < oo. Function G {x) is called tempered, if there 
exists some constant a such that 

\G {x)\ < \x\'' \ x-.±oo. (1) 

Fourier transform of any tempered distribution is a tempered distribution as 
wellHISlEI. 

The standard procedure of transition from Fourier series to Fourier integral 
is based on infinite increasing of r (see e.g. [3IH]). Usually one starts from 
locally summable function F {x) and constructs a truncated function F {x) as 

^^""^"I 0; if |x|>7rr 

Then F (x) is extended periodically (with the period 2ttt) on the whole real 
axis X, i.e. F (x) — > F (x/t) where^ 

F{^)^F{^ + 2n) (3) 

and _ 

F{x/t)^F{x)^F{x)- (4) 

for —TTT < X < TTT. 

This procedure may be easily generalized for a case when F {x) , F {x) , and 
F{x/t) are distributions. It is clear, that since generally distributions are not 
defined pointwise, the periodicity condition Q should be interpreted as. 



Fcj,^) = {F<j,^+2^) (5) 

where 0^ is arbitrary test function 0^ G with support located in the interval 
((^ — e, + e) with arbitrary small, but finite e. 

It is evident that the periodic distribution F{x/t) may be obtained as the 
result of the cyclization procedure ^| defined as 

oo 

F{x/t)^Y:F{x)^ ^ F{x + 2TTmT) (6) 

m— — oo 

Periodic distribution may be represented as the Fourier series 

oo 

F {x/t) = ^ F„ exp {mx/T} (7) 

n— — cxD 

where Fourier coefficients Fn are given by the standard expression 

1 ~ 

Fn — / F (s/r) exp {— mx/r} dx (8) 

^ Below the tilde always marks periodic functions with the period 2tt. 
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Fourier coefficients can be also computed as 

1 f'^ — 

Fn — / F (a:) exp {—inx/r} da;, (9) 

Although correspondence between Fn and F{x/t) is biunique, there are no 
biunique correspondence between F„ and F {x) in a general case (see Appendix 

I)- 

Lastly, we wish to stress that all mentioned error and ambiguity problems 
are of no particular interest, if the considered transition is used only to define 
Fourier integral, as it is e.g. in |7||H1. When such definition is done, even for some 
restricted class of functions, further extension of Fourier transform on a wide 
class of functions and even on distributions needs no reference to the original 
Fourier series. Bijection is needed only for the function and its transform, given 
by the Fourier integral. However, there are cases when the transition from 
Fourier series to Fourier integral is of particular interest, e.g. when r — > oo 
corresponds to the lifting of regularization. Hence the degree of ambiguity in 
such transition deserves a more detailed consideration. 



2 Transition from Fourier series to Fourier inte- 
gral 

After [H] we split the region of summation — oo < n < oo into 'bursts', i.e. 
intervals rfc — T5k/2 < n < rk + TSk/2 and approximately^ replace the sum by 
its average value in each interval 

^ Tk+T5k/2 

— ^ exp {inx/r} Fn ^ t [exp {inx/r} Fn]n=rk = exp {ikx} tFut 

n—Tk — T5k/2 

(10) 

Some details of partition into 'bursts' may be clarified, if we note that after 
formal change n — > fcr we get from 

1 r'' ~ 

— / F {x/T)eyip{—ikx} dx — fkT6k (11) 

271" J-nr 

with definition 

/fe = TFkr (12) 

It is assumed in jH|, that when proceeding to limit t ^ oo, the expression 
was converted into the inverse Fourier transform 

1 r'' ~ 

— lim / F {x/T)exp{—ikx} dx — fk (13) 

2tt t^ocJ_^^ 



■^Although the term 'approximately' isn't specified in 0, one may interpret it so that the 
approximate equahty is assumed to become an exact one with the vanishing Sk and l/r. 
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From we conclude, that it may happen only if for large enough r, we impose 
the condition 

Tdk = 1 (14) 

it means that 'bursts' are defined as rk — 1/2 < n < r/s+1/2, so that each 'burst' 
contains only the sole term exp {ikx} Fkr , where fcr is still an integer number. 
Relation H12|) is nothing but a new definition so it gives nothing new concerning 
the conditions under which right-hand member in (|1U|) may be treated as an 
integrand of the Fourier integral in variable k. 

For lack of something better we suggest to claim the fulfillment of approxi- 
mate relation'^ 



Fn exp {in(p} ~ / Ft exp {itip} dt; < e < 1 

Jn+e~l 



(15) 



which directly provides the desired resuh 



n— — oo 



Fn exp {in(p} ^ I Ft exp {itif} dt = tF^t exp {ikx} dk (16) 

_ ^ , J —oo J —oo 

/oo 
fk exp {ikx} dk (17) 
-OO 



If we define 

/oo 
fk exp {ixk} dk (18) 
-oo 

then ()17|l may be rewritten as 

F{x/t)^F{x) (19) 
With the definition (|12|l the approximate relation H15|) acquires the form 

fk exp {ikx} ^ jr ft exp {ixt} dt (20) 

"1^ Jk-(l-e)Sk 

Relation ()20|l resembles the result of the mean value theorem, however ft should 
not obligatory obey the restrictive conditions of such theorem because (pn|l must 
be satisfied only for infinitesimal Sk and only approximately. 

Making allowance for ^ , the inverse Fourier transform may be written with 
the same accuracy as 

1 

fk = Tr ^ (^) {-ikx} dx (21) 
In other words, ordinary transition may be done by the formal substitution 

n^kr; ^ ^ j ^" = ^-^'^^ ^^2) 



^The error of such approximation will be estimated later. 
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In this paper we try to clarify under what conditions substitution H22I) leads to 

lim F{x/t) ^ F{x) (23) 

r — *oo 

which we treat as an exact result. 

It is clear that nontrivial result appears after proceeding to limit r — > cxo in 
H12|) and (|23|l . only if Ft and F (x/t) have specific properties. Some of them are 
discussed in Appendix II. 



3 Ambiguity in choosing of extension 

In the distribution theory Fourier series exist even for such distributions that 
obey neither Riemann-Lebesgue lemma nor the Dirichlet condition. As it is 
proved in [B', Fn exist and are called Fourier coefficients, if F {(p) is locally 
summable (i.e. at any finite interval). Distribution theory allows to sum series 
in which F„ is a tempered function of discrete argument n. Therefore, an appro- 
priate tempered distribution may be taken, as an extension Fn for noninteger 
n. For instance, if Fourier coefficients are given by 

Fn = ^ Q. if ^j<0 ; X = RcX = const (24) 
straightforward extension n ^ t gives the distribution 

If extension is allowable, then is allowable too under the condition 
that Fj^' — Ff'^ = LOt turns into zero for all integer t — n, i.e. ujt & jlOj . 
However, when a class of the sought distribution is specified, the choice of 
extension may be essentially restricted. If, for example, it is required that the 
target distribution must have continuous derivative up to m-th order in some 
specified area (e.g. in the vicinity t — 0), then straightforward extension (|25|l is 
unsuitable if m > A. 

There exists an extension of Fn with exceptional analytical properties, namely 
the entire function 

1 r'' ~ 

Ft^ / F{x/T)exp{-itx/T}dx (26) 

which may be obtained from expression IHll. When according to problem sit- 
uation the extension H26|l is prescribed, then even if Ft breaks the conditions 
of Carlson theorem (see e.g.jHI), it may be shown that such extension JO] is 
unique. 

On the other hand, if F {x/t) is known, F [x) may be simply computed with 
H23I) . In many interesting and more common situations technical difficulties 
prevent summarizing series in (|7|l. 
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Formally, an extension may be reconstructed from any other, without 
direct computation of series in (|7|). Indeed, the substitution of Q in (I26|l gives'*: 



Ft^— exp{-it(p} ^ Fn exp {irup} dif ^ ^ F„ 

' ' — — rv~, n — — rv~\ 



sin7r(n-t) ^^^^ 



(n-t) 



where -Fi is an entire function. Unfortunately, the formal expression (|27|l adds 
nothing to computation of fk or F (x), since (|22|l with (^I] 5(58)) gives 

-,00 • i 1 n\ 00 

1 sm TT fc — - r .r-^ / 71 \ , , 

A-^:^ E ^« E (28) 



Tl— — 00 n— — 00 

and Fourier transform of (|28|l brings us back to 10. 

If we take an exact summation of series in H27f) we, undoubtedly, obtain a 
practicable expression. Unfortunately, as a rule the sum computation in ()27|l is 
as difficult as the summation of series in Q , so we consider the simple example 

i sinh i/TT cos 7rn . 
TT n + iiy 

where v is positive noninteger number. From (|27|l we find that the corresponding 
entire function is 

^ siny(t + z.) 

As it should be, the difference of Ft obtained from H29|) by the straightforward 
extension n ^ t and Ft from (|3U|) 

cosh :/7r sin iTT ^ 
TT r + zz^ 

For sure, Fourier series either with the coefhcients given in (|29() or those given in 
H30I) converge to the same periodic function exp{z/^}, which is equal exp{i^(^} 
in the interval — tt < (/j < tt. 



4 Error and ambiguity introduced by formal tran- 
sition 

Although procedure proposed in [71 |S] is entirely appropriate to define the 
Fourier integral, being applied to transition from Fourier series to Fourier in- 
tegral, for a wider class of functions it introduces appreciable difference be- 
tween F (x) obtained after formal transition (|22|l and an exact one: F (x) = 
limT-^oo F (x/t). To estimate the error we apply the Poisson formula 
00 00 -.00 

J2 <i>„exp{*n^}= / e'('^-2™)*^^^^ (32) 

n—~oo m— — oo 

*The Abel-Poisson regularization is implied, so the change of the order of summation and 
integration is legal. 
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to the right-hand member in ||7)|. It gives 

~ °° f°° fx 1 

F{x/t)— / FnCxplin 2TTinin> dn (33) 

that with (|12|l may be rewritten as 

— roc 

F(x/t)^ / fkexp{ik{x + 2TTnT)}dk (34) 

and, taking into account (fTH|l . we finaUy get 

oo 

F{x/t)^ F{x-2TrnT) (35) 



n— — oo 



Since both (|12|l and (|18() are simply definitions and include no approximations, 
1)35(1 is an exact expression. Comparing it with (|19|l we see that F {x/t) differs 
from F {x) in (QHl) by 

A (x) = ^ (x - 27r™) (36) 

With T — > cxj correction A {x) must disappear for any physical value F {x) . It 
is not true, however, for auxiliary values, e.g. for gauge fields. 
For instance, from 

Ft=LUt = {1- exp {2iTrt}) Rt e n (37) 
which turns into zero^ for all integer t, we obtain after the formal substitution 

/fc = (1 - exp {i2nkT}) rk (38) 

where 

Tk = TRkr (39) 

so Fourier transform of (|38|) gives 

F{x) =r (x) +r{x + 27rfcr) (40) 

hence if 

/>oo 

r{x)= / rke'^'^dk (41) 



decreases with x — > oo, then for r — > oo and any finite x we get (x) = r (x), 
despite the fact that corresponding Fourier series is identically zero. 

Nonetheless, in a case when extension (I26|l is prescribed, no ambiguity arises. 
Indeed, if we apply the formal substitution to extension (jJEl, we get for r oo 

fk =TFkr^— F {x/t) e-'^^dx (42) 
27r . 1^^ 



^The distribution Rt is assumed not too singular at integer t to spoil that condition. 
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On the other hand, from the definition (|13|1 we obtain 

p7rT-\-2n7iT 

27r 



E ^ / F{x)e-^^^dx (43) 



n— — cxD 



2 

that with (|42|l evidently leads to the exact result (|35|l . 

5 Asymptotic series expansion 

As it can be seen from (I34f) . one may estimate the correction A (x) by com- 
puting the limit of exp {iA;27mT} for r — > oo. In |12l I13| the method of 
asymptotic series expansion of exp {i/cr} e*'^'^/fe for r — > ±oo was suggested. 
In particular, it is shown that e'^'^/fc vanish in such limit, if fk is analytical 
in the vicinity of A; = 0. In our case it means the vanishing of A (a;). An 
exhaustive investigation of the power type distributions, i.e. distributions of 
a form / {x) ~ (x ± iO)'^ In'" {x ± iO) has been undertaken in [131 111) . In [TH 
this method was extended over a family of distributions and the case when the 
regularization term in f^. depends on t was considered. It was shown that 



n=0 

where 



/•oo /"O 

/^+) = / i^(i)exp{itA:}di; fjr''^- F (t) exp {itk} dt (45) 

and 

27r ''-^ n — m to 

m=o ^ > 

where (t) is n-th derivative of the Fourier transform of 
In particular, for A = — 1 we get 

e-*^-" (fc + i^)"' = Tie^"^ E,T=o if r ^ ±oo (47) 

and ^ 

^fc + i-) = if T^icx) (48) 



so returning to the example considered above, we see that after formal substi- 
tution (|^ from (j^ we get for r —> oo 

= l£i!!iii^£^^ _ rF,. = /(^) . sinh (y.) e~-5 (k) + O (1/r) (49) 
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and from we obtain 

cosh z^TT sin nTT (^) /s) , . x -7r,^xn\ i ^ /i / \ f^n\ 

' = — ^j-j— ^ ^fer ^ fk - cosh (z^Tr) e 5 (fc) + O (1/r) (50) 

so 

A = /i^^+/i^^='5(fc) (51) 

and, consequently, 

F{x) = l (52) 
that, indeed, coincides with the exact result 

F (x) = lim F (x/t) = lim exp < ii/x/r f = 1- (53) 

6 Conclusions 

The transition from Fourier series to Fourier integrals is considered and errors 
introduced by the formal substitution are estimated. Also the ambiguity 
caused by transition from discrete function to continuous one is considered. 
Conditions under which such ambiguity disappears are suggested. 

7 Appendix I. 

As it is shown in 10 , the distribution 

Af^(x) = p(^-7r) Gf!' (54) 

have only trivial periodical extension 

AF(a;/r) ^E(AF(a;)) =0, (55) 

so both F (x) and Fp {x) = F (x) + AF {x) after the cyclization procedure I© 
give the same F {x/t). It also means that inverse procedure F (x) = {x/t) 
is ambiguous, i.e. F (x) cannot be reconstructed uniquely for a given F{x/t). 
Since F p [x) should vanish outside —tit < x < ttt, distributions p ± tt) may 
differ from zero only at the points x = I^ttt, i.e. p (^ ± tt) have point support. 
As it is known (see e.g. 01 E]), distribution with point support is finite linear 
combinations of Dirac (5-functions and its derivatives, so one may choose 

M 

AF{x) = ['^^"^ (7 " " '^^"^ (7 + ' 

m— 

with arbitrary finite constants Cm and M . 

As it can be seen from (I55|) no contribution may come from AF {x) to Fourier 
coefficients, in other words 



AFt = — / AF (x) e-"^/^da; 



(57) 
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must turn into zero for all integer t ~ n. To guarantee this, we must remove 
some uncertainty in integration J^™^ ± tt) on [—ttt,ttt]. Recall that Dirac 
(5- function is defined as the functional/^ (y) (5 (y) dy — 0(0), for any test 
function (j) (y) continuous in the vicinity y ~ 0. However, there exists no such 
4> (y) with support located in the interval < y < oo, i.e. 

/CO /-oo 
<l>{y)S{y)dy= / ^{y)Siy)dy (58) 
oo Jo 

with <j) (0) 7^ 0. Certainly {(p, 6)^0 does not look winning at all. Even if we 
put (2/0 = then 

_ 1 *^ 

AFt = Ihn [e'^'e (e) - e-'^*0 (-e)] — ^ (zt)'" c„ (59) 

m=0 

SO we only shifted the trouble to 9 {x), which is undefined at a; = 0. To be more 
exact, value 9 (0) may be chosen arbitrarily. 

Indeed, let as consider some function D (x) which possesses arbitrary values 
on a countable set of x and is equal to zero on a complement of such set. We 
assume also that D (x) does not include linear combinations of Dirac (5-functions 
and its derivatives. In this case the functional {D, 0) = for any of the test 
functions (f). The distribution theory expressly makes no distinction between 
distributions, which differ by D(x), if D (x) satisfies the above condition. In 
other words, on a countable set of points one may arbitrarily choose values of 
any distribution, if at these points such distribution can not be presented as 
Dirac (5-functions and its derivatives. 

In particular, one may choose so called 'symmetric' definition 

OW^l.^n C'"-^';""'--' ). (60) 

that gives 9 {0) — 1/2. The 'symmetric' definition is especially convenient, if 
distribution is represented as the Fourier series, because, according to Fejer 
theorem (see e.g. [15J), if Fourier series is convergent at </? + and ip — 0, it 

converges to F (ip) — (^F [ip + 0) + F {ip ~ 0)^ /2. The 'symmetric' definition 

allows to get for AFt an acceptable result. 

. . . M 

AFt^z^^ V(zt)™c„ (61) 

TT — ' 

Nevertheless, that or another choice is merely a matter of convenience, because 
this problem lies outside the scope of the distribution theory. For example, in 
^J8. 1(732) another definition is taken 

/ e-P'^S (x) dx^ S{x)dx = l (62) 
Jo Jo 

i.e. 9(0) = 1. 
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8 Appendix II. Some remarks about analytical 
properties of Ft and F{x/r) 

As it already mentioned, after proceeding to limit r ^ oo the nontrivial result 
may appear, only if Ft and F {x/t) have specific properties. Indeed, if F (x/t) 
depends on x and r only through x/t, then limT-^oo F {x/t) — F (0) for any x 
where such limit exists, unless x is the infinite point. In particular for any finite 
Fourier series, one may write 

F {x/t) = F„ exp {inx/T} ^ F„ = F (0) = const (63) 

n=-N_ n=-N_ 

for any fixed N±. 

In proceeding to limit t — > oo the positions of singularity of TFkr and F {x/t) 
are essentially rearranged. For instance, if all singularities of Ft are located in 
the finite area \t — to\ < R, where to and R may be arbitrarily large, but finite, 
then for T ^ oo all corresponding singularities of fk are moved to the infinitely 
small vicinity of fc = 0. Indeed, in our case Ft may be expanded in Laurent 
series 

for any t which obeys \t — to\ > R. Therefore, after the formal substitution (|22|l 
we get 

Ft^fk^ lim TFkr = lim V j-, , , (65) 

T^oo T-too — ' (k — tn T) 

Ai 



k ~ ie signum (Im to) 



Less trivial and more realistic result may be obtained only when r-dependence 
of A„j is allowed. We exclude the case when Am increase faster than r™^^ with 
T — > oo, because it makes fk infinite for all k. Therefore, either A„i — > CmT™'~^ 
or corresponding terms vanish. Thereby, we may write 



oo 



fk = lim " m = 1™ 7 Wm- (66) 

T^oo {k — tn t) ■^-^ (fc — z£ signum (Imtn)) 

m—1 ^ ' ^ m—1 ^ o \ y/ 

So if series (|66|l converges, fk became an analytic function in a whole complex 
plane fc, except the point fc = 0. 

Let now F {ip) be regular in the area \ip\ < c, where c is arbitrary small, 
but finite, then F {x/t), and consequently F {x), may have singularity only at 
\x\ > CT, so after r oo distribution F{x/t) converts into entire function. 
When F {ip) depends not only on (p, but on r as well, i.e. F = F {(p,T), and if 
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F (t/?, r) is analytical in a ring r_ < \ip\ < r+, i.e. 

oo 



n— — oo 



with r± — lim„^±oo _F„ (0, r) , such {(p, t) may have nontrivial behavior 
for infinitely increasing period under the conditions which may be easily found. 
Indeed, within such ring F (ip, r) may be presented as Laurent series and for a 
reason similar to the one mentioned above we obtain for r — > cxd 

K (0, r) ^ t"$„; = const (68) 

so the distribution 

oo 

<^>{x) EE \im F (-,t) ^ y $„x-"-i (69) 

T^OO \T / ^ — ^ 
n— — oo 

is analytic in a ring i?_ < \x\ < i?+ with R± = lim„^±oo |*i>n| " ■ 
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